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Abstract 



O 

I The DELPHI Collaboration has recently reported the measurement of J/ ip production in photon- 

photon collisions at LEP II. These newly available data provide an additional proof of the impor- 
tance of colored cc pairs for the production of charmonium because these data can only be explained 
by considering resolved photon processes. We show here that the inclusion of color octet contribu- 
tions to the J/ip production in the framework of the color evaporation model is able to reproduce 
this data. In particular, the transverse-momentum distribution of the J/ip mesons is well described 
by this model. 



PACS numbers: 13.60.Le, 14.40.Gx 



*Elcctronic address: eboli@fma.if.usp.br 
'''Electronic address: gregores@ift.unesp.br 
^Electronic address: mizuka@fma.if.usp.br 

Typeset by REVTeX 



I. INTRODUCTION 



DELPHI Collaboration recently released preliminary measurements of the transverse mo- 
mentum spectrum of J/ip mesons produced in 77 collisions at LEP . This new data 
allow further tests of models for charmonium production. We show here that the Color Evap- 
oration Model (CEM) reproduces these new results using the same single non-perturbative 
parameter that has been obtained from previous analysis of charmonium photo- and hadro- 
production. These newly available data provide an additional proof of the importance of 
colored cc pairs for the production of charmonium, as the data on this region can only 
be explained by considering resolved photon processes, which forms colored cc pairs in the 
leading order. The CEM for charmonium production incorporates these colored pairs into 
the total yield of charmonium in a very simple and economical way. 

The Tevatron data 0, 0] on charmonium production at high changed the way we 
understand charmonium production. The presently successful models are based in two key 
considerations: i) onium production is a two-step process where a heavy quark pair is 
produced first, followed by the non-perturbative formation of the asymptotic states, and 
ii) color octet as well as singlet cc states contribute to the production of charmonia. These 
'eatures are incorporated in the Non-Relativistic QCD (NRQCD) factorization approach 
^,EJ, in the Color Evaporation Model [^HsHil, and in the Soft Color Interaction Model jloj|. 

The Color Evaporation Model simply states that charmonium production is described 
by the same dynamics as DD production, i.e., by the formation of a cc pair in any color 
configuration. Rather than imposing that the cc pair is in a color singlet state in the short 
distance perturbative processes, it is argued that the appearance of color singlet asymptotic 
states solely depends on the outcome of non-perturbative large distance fluctuations of 
quarks and gluons. These large distance fluctuations are considered to be complex enough 
for the occupation of different color states to approximately respect statistical counting. In 
fact, it is indeed hard to imagine that a color singlet state formed at a range m^ 1 would 
survive to form a ^ at a range Aqq D . Although far more restrictive than other proposals, 
CEM successfully accommodates all features of charmonium production ^ , 12 , k| . 

The CEM predicts that the sum of the production cross sections of all onium and open 
charm states is described by 

1 f 2m o da c5 
Oonium = r / dM c5 — — , (1) 

y J2m c aMcc 



and 

r2m D dorr f , , , da, 



a opcn = - / dM c5 — f- + / dM c5 — f- , (2) 

9 J2m c dM C c J2m D dM c5 

where M c5 is the invariant mass of the cc pair. The factor 1/9 stands for the probability 
that a pair of charm quarks formed at a typical time scale l/M^ ends up as a color singlet 
state after exchanging an uncountable number of soft gluons with the reaction remnants; 
for further details see Q. One attractive feature of this model is the relation between the 
production of charmonium and open charm, which allows us to use the open charm data 
to normalize the perturbative QCD calculation, and consequently to constrain the CEM 
predictions. 

Up to this point, the model has no free parameter in addition to the usual QCD ones. 
In order to predict the production rate of a particular charmonium state, let us say a J/ip 
meson, we must also know the fraction p^ of produced onium states that materialize as this 
state (J/ip), 

°V> = Pip ionium • (3) 



In its simplest version, the CEM assumes that p^ is en ergy and process independent, which 

is in agreement with the low energy measurements Notice that p^ is the solely free 

parameter of the CEM, making this a very restrictive framework. From the charmonium 

ri 

photo-production, we determined that p,/, = 0.43-0.5 |8|, a value that can be accounted for 
by statistical counting of final states [lOj. The fact that all ip production data are described 
in terms of this single parameter, fixed by J Jib photo-production, leads to parameter free 
predictions for Z-boson deca y ra te into ip [la] , and to charmonium production cross section 
at Tevatron and HERA jisl ll^. as well as in neutrino initiated reactions [2^. 



II. RESULTS 

The differential cross section for the inclusive process e + e~ — > e + e~77 — > J/ipX is 
= E JfJJ d y +d y~ dx Adx B f 1 / e +(y + )f l/e -(y-)F Ah (x A )F Bh (x B ) d ^ , 

(4) 

where / 7 / e ± is the bremsstrahlung photon distribution from an electron/positron. We de- 
noted the parton distribution function of the photon by Fa{b]/~/(xaib}), where xa[b] is the 
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fraction of the photon momentum carried by the parton For direct photon inter- 

actions (A[B] = 7), we have Fa\b]/*/{%A[b]) — $( X A[B] — !)• We considered an average 
electron-positron center-of-mass energy 2E e = 197 GeV. We also applied the experimental 
J/i/j rapidity cut —2 < r)^ < 2, and imposed that the 77 center-of-mass energy satisfies 
W J7 < 35 GeV, where W 77 = 2E e ^y~^r . 

In our calculation, we employed the Weizacker- Williams approximation for the photon 
distribution 



/ 7 /e±(y) 



2ti 



1 + (1 ^ )2 log(g-W/ 1 



(5) 



v \ O 2 I e \ O 2 O 2 

t> \ mm / \ max ^ mm , 

with Qmin = m ly 2 /( l ~ V), and Qmax = (E e 6) 2 (l - y) + Q 2 min . Here, the fraction of the 
parent e ± energy (E e ) carried by the photons is y (= E y /E e ), and 6 is the angular cut that 
guarantees that the photons are real. We used 6 = 0.032 radians, as determined by the 
experiment. 

The inclusive subprocess cross section a(AB — > ipY) was calculated using the CEM; see 
Eqs. (JTJ) and ©• The partonic subprocesses contributing to J/ip production are depicted 
in the Table HI Notice that both direct and resolved photons contribute to charmonium 
production in the CEM. We evaluated numerically the tree level helicity amplitudes of the 
subprocesses displayed in Table HI using MADGRAPH Q and HELAS [22] packages. The 
adaptative Monte Carlo program VEGAS H was employed to perform the phase space 
integration. 

In the framework of the CEM, the evaluation of the photon-photon production cross 
section contains only the free parameters appearing in the perturbative QCD calculation of 
the subprocesses presented in Table HI since the CEM free parameter can be fixed at the 
value extracted from the photo-production of J/i/j, i.e. = 0.5 [sj]. We used the leading 



Direct Once Resolved Twice Resolved 

77 ccg ~fq(q) ccq(q) qq -> ccg 

19 -> ccg gq(q) ccq{q) 

99 -> ccg 

TABLE I: Subprocesses contributing to J/tp production in 77 collisions. Here q stands for the 
light quark flavors u,d,s. 
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FIG. 1: Uncertainty on the p\ differential cross section originated from different choices of the 



renormalization scale /jr. In (a) we chose /jlr = £m c while in (b) /ir = ^Jm^ + Pj<- The shaded 



band was obtained by varying \ < £ < 2. We fixed 



m. 



1.3 GeV, and used the GRS-G parton 



density function in both figures. 

order GRV-G Q and GRS-G |2j parton density functions as provided by CERN PDFLIB 
package with the partonic subprocess center-of-mass energy as factorization scale /j,f — V§- 
We verified that our predictions do not vary significantly for other choices of the factorization 
scales, e.g. /if = ~ \/I and fip = 2 y/S. We also verified that the results are very similar for 
the GRV-G and GRS-G parton distributions (see Table ITT]) . The strong coupling constant 
was evolved in leading order considering four active flavors and A.q CD = 300 MeV, while the 
charm quark mass was varied between 1.2 and 1.4 GeV. 

In order to access the theoretical uncertainties in the lowest order CEM calculations, 
we analyzed the predicted J/ip transverse momentum spectrum for different choices of the 
renormalization scale (hr) ■ We present in Fig. the predicted spectrum obtained for 
ur = £m c with | < £ < 2 and m c = 1.3 GeV, as well as the DELPHI experimental results 
|l|, |2q . We can see from this figure that CEM describes well the shape of the distribution, 
despite the large uncertainty in the absolute value of the differential cross section. Notice 
that we are only changing a global factor (as) for this choice of fiR when we vary £. Figure 



HJd displays the p\ spectrum for [i R = £ym^ + p\ with \ < £ < 2 and m c = 1.3 GeV. For 
this choice of fiR the uncertainties in the the p\ distribution are smaller than for the previous 
choice of hr. However, the shape of the p\ spectrum changes and the CEM prediction seems 
to diminish faster at large p\ than the data. 
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FIG. 2: Differential cross section as function of the squared transverse momentum of the J/ip. The 
shaded band shows the theoretical prediction obtained by varying the charm mass (m c = 1.3 ± 0.1 
GeV). We explicitly show the contributions from direct, once resolved and twice resolved cross 
sections for m c = 1.3 GeV. 

In Figure 121 we display the contributions to the J/ip Vr spectrum arising from direct, once 
resolved, and twice resolved processes. These distributions were obtained using the GRS-G 
photon parton densities, /i^ = m c and m c = 1.3 GeV. As we can see, the once resolved 
processes are responsible for the majority of the events (~ 85%) while direct and twice 
resolved processes account for less than 15% of the total cross section. The most important 
process is jg — > ccg. We also present in this figure the uncertainties associated to the charm 
quark mass; the shaded band represents the sum of all contributions taking m c = 1.3 ± 0.1 
GeV. Notice that the largest uncertainties in the CEM prediction originates from the choice 
of the renormalization scale. This is quite expected since we are performing our calculation 
in lowest-order perturbative QCD. We summarize our results for the total cross section in 
Table HH 

In order to further _ ou r resets with the recent* puhhshed DELPHI resnUs QQ. 
we evaluated the dependence of the total J/ip yield on the minimum transverse momentum 
for = 197 GeV. The result is presented in Fig. As can be seen from this figure, the 
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Minimum pf (GeV 2 ) Energy (GeV) 



FIG. 3: Total cross section as function of the minimum squared transverse momentum (a) and the 
e + e~ center-of-mass energy (b). In (a) the solid line stands for the DELPHI measured total cross 
section while the dotted lines indicate the experimental error of this quantity. We varied the charm 
quark mass as m c = 1.3 ± 0.1 GeV to estimate the theoretical uncertainties. We used \fs = 197 
GeV for the minimum transverse momentum dependence (a) and imposed p 2 - > 0.25 GeV 2 for the 
center-of-mass energy dependence (b). The remaining parameters are the same as for Figure |2J 

choice of QCD parameters we used in this analysis provides a very good description of the 
existing data, reinforcing our confidence on the predictive power of of the color evaporation 
model. Figure displays the CEM predictions for the J / ip production cross section as a 
function of the e + e~ center-of-mass energy (y/s). Here we assumed that ip\ > 0.25 GeV 2 , 
ix r = m c with m c = 1.3 ±0.1 GeV, and we used the GRS-G set of parton distribution 
functions. As expected, the total cross section grows with the center-of-mass energy due to 
the increase in the photon-photon luminosity. We verified that contributions of direct, once 
resolved and twice resolved processes are in the same proportion of the results presented for 
\fs = 197 GeV; see Fig. EJ Taking into account the planned luminosity of the future e + e~ 
colliders, we can easily foresee that it will be possible to extract very precise data on the 
photon-photon charmonium production in these machines. 

III. CONCLUSION 

In this paper we showed that the Color Evaporation Model for quarkonium production 
correctly describes DELPHI data on J/ip via photon-photon collisions. Due to the rather 
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Parameters Cross Sections (pb) 
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L/ J.X V V Ij 


Oticp RpsoI vpH 


Twirp T?PF;olvpf] 

J- VV J. W J. LvOWX V WL 
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ORV-G 


1 3 
x .o 


1 n 

x . \J 
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1 79 


1 3 3 

XO. (J 


1 00 
x . uu 


16 1 

1U.1 
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1 3 


1 n 


n 


1 72 


1 3 


94 


15 7 

XfJ . 1 


GRS-G 

v.l 1 v.l 


1.2 


1 n 


n 


2 75 


21.8 


1 73 


26 3 


GRS-G 

v.l 1 V,kJ v.l 


1.4 




n 


1.02 


7 5 


51 


9 07 


GRS-G 


1.3 


0.5 





3.26 


46.8 


6.42 


56.4 


GRS-G 


1.3 


0.5 


1 


2.42 


28.4 


3.15 


34.0 


GRS-G 


1.3 


1.0 


1 


1.44 


9.65 


0.61 


11.7 


GRS-G 


1.3 


2.0 





1.17 


5.99 


0.29 


7.45 


GRS-G 


1.3 


2.0 


1 


1.03 


4.84 


0.22 


6.08 



TABLE II: Cross sections for direct, once resolved, and twice resolved production processes for 
v\ > 0.25 GeV 2 using different sets of parton distribution functions, charm masses, and renormal- 



ization scales /j_r = £,ym^. + 0Pj>- 

large uncertainties in the data, its is not possible to use them to discriminate between the 
different proposed mechanisms for charmonium production. As far as the DELPHI data are 
considered, the NRQCD 2(| and CEM frameworks present equivalent results. 

Considering that the CEM is also successful in describing the photo- and hadro- 
production of charmonium, we conclude that this model gives a robust and simple pa- 
rameterization of all charmonium physics. Moreover, 77 reactions provide a clear proof of 
the importance of colored cc pairs to the production of charmonium, since the data on this 
reaction can only be explained considering resolved photon processes, which lead to colored 
cc pairs. 
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